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We develop a diagrammatic perturbation treatment to calculate the zero-temperature equation 
of state of the dilute gas mixture of a single spin component Bose-Einstein condensate (BEC) and 
QO | a normal Fermi gas of indistinguishable (single spin) fermion particles. We find that the mean-field 

description breaks down near the mechanical instability related to the phase separation phenomenon. 
Our analysis shows that the instability is caused by the competition of the usual short-range and 
fermion-mediated boson-boson interactions, which result in a boson compressibility that diverges. 
In the low BEC-density limit, we show that the diagrammatic analysis simplifies, we sum part of 
the higher order diagrams, and we discuss the effects of other higher-order contributions. 
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I. INTRODUCTION 



Cold atom technology provides a novel laboratory for the many-body study of quantum liquids. In this paper, we 
focus on the low-temperature cold atom mixtures of dilute boson and fermion quantum liquids. The simultaneous 
trapping and cooling of fermionic and bosonic atoms has become a routine ingredient of the experimental cold atom 
repertoire. The original motivation for trapping these quantum gas mixtures was sympathetic cooling [l], 0, In the 
& ■ quantum degenerate regime, the standard evaporative cooling technique-the removal of the most energetic particles-is 
more efficient with bosonic than with fermionic atoms. Evaporatively cooled bosons, often a Bose-Einstein condensate 
(_J . (BEC), can subsequently cool the fermion particles by thermal contact. 

In addition to its sympathetic cooling use, theorists also have pointed out intriguing prospects of cold atom technol- 
ogy for many-body studies: As the elementary BEC excitations are acoustic phonon modes, fermion-BEC mixtures 
consist of fermions interacting with a phonon fluid, thereby representing a large class of systems. A cold atom feature 
that is highly unusual in this broad context is the continuous tunability of the fermion-boson (and hence fermion- 
phonon) interactions, obtained by adjusting the strength of a homogeneous, external magnetic field-the Feshbach 
| resonance. In general, the unprecedented accessibility, novel knobs, and probes suggest that in their role of proto- 
types of fcrmion-phonon fluid, the cold atom fermion-BEC mixtures can be used to break new ground by exploring 
polaron-like self-localization [4] or by measuring the rate of heat-exchange between strongly coupled fermion-phonon 
. fluids (temperature relaxation). 

From the perspective of quantum liquid physics, a particularly relevant phenomenon is the phase-separation transi- 
tion that the fermion-BEC mixture is expected to undergo as the fermion-boson interaction strength is increased [H . 
This transition is the analogue of the phase separation of condensed 3 He — 4 He fluid mixtures studied in traditional 
low-temperature physics Q. As we describe in Sec. [Til we expect the cold atom phase-separation transition to be 
\Q • first order Q, as was observed in the condensed-matter helium mixtures. Accessing this transition by varying the 
densities and/or the interaction strengths would allow the unusual, though highly relevant, experimental exploration 
of first-order (zero temperature) quantum phase transitions. As a first-order transition, its dynamics can involve 
nucleation or a spinodal decomposition. At sufficiently low temperatures, the nucleation process cannot be triggered 
^ , by thermal activation and has to proceed via a many-body tunneling process. Hence, the cold atom fermion-BEC 
mixtures promise the prospect of many-body tunneling studies, as well as investigations of first-order quantum phase 
rN | transitions [Ic| |. 

The focus of our work in this paper is the role of mediated interactions-the modification of the interactions of 
particles in one fluid caused by the presence of the other fluid. This phenomenon is of fundamental importance: 
fermionic particles interact via boson-field mediated interactions in the standard model. In the phase separation of 
low-temperature condensed helium mixtures, mediated interactions play a crucial role [Xlj | . In those systems, the 
condensed 3 He fermion liquid remains "normal" (i.e., does not become superfluid) in the experimentally accessible 
regime, whereas the 4 He boson fluid can take on a BEC-like superfluid state. Boson mediated interactions are expected 
to make the fermions superfluid |l2] as well, but at temperatures that are an order of magnitude lower than can be 
reached with present cooling technology. The condensed 3 He — 4 He mixtures undergo phase separation when the 
fraction of fermion particles exceeds a minimal value: the system then breaks translational symmetry by forming 
local regions in which only fermion particles reside. While our understanding of theses systems remains limited by 
the difficulties involved in accounting for strong interaction effects, theorists understood early on that the phase- 
separation phenomenon is triggered by mediated interactions [111 ]. Mean-field studies of the cold atom fermion-BEC 
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mixtures indicated that these gas mixtures would similarly undergo phase separation under experimentally attainable 
conditions. Hence, it should be possible to study mediated interactions in a cold atom environment that is amenable 
to a first-principle description. 

Mediated interactions stem from a particular type of correlations: the local density of the second fluid is altered 
near the position of a particle "1" of the first fluid, thereby changing the potential energy experienced by another, 
nearby particle "2" of the first fluid. However, most cold atom phase separation studies have been carried out in the 
mean- field approximation, which neglects all correlations. Other studies neglected dominant contributions [15, [Til. [l5| . 
The zero-temperature equation of state and, hence, the density profile calculated within the mean-field equation of 
state cannot be expected to give a quantitatively correct description. 

In Sec. we derive the equation of the state of the above system within the mean-field approximation. In Sec. IIIIl 
we describe the fermion-mediated boson-boson interactions starting from linear-response theory. Then, we develop a 
general perturbation formalism that accounts for the relevant correlations in the following four sections, Secs. lIVHVIII 
Our second-order calculation indicates that the instability associated with the spinodal decomposition is not a saddle 
point instability, as predicted by the mean-field model, but is caused by the divergence of the BEC compressibility. 
In the low-density BEC limit, we find that the perturbation treatment simplifies and we show that higher-order terms 
significantly lower the fermion density of the instability point and make the instability line dependent upon the BEC 
density. From a partial summation of diagrams, we find an analytical expression for the equation of state and estimate 
the effect of other contributions. 



II. COLD ATOM BACKGROUND AND MEAN-FIELD DESCRIPTION 



In this section, we describe the ground state of a dilute gas mixture of single-component fcrmionic atoms and single 
spin projection bosonic atoms. These fcrmions (bosons), interacting via short-range atom-atom interaction potentials, 
are indistinguishable particles with mass tuf (ms), distributed in space with average particle density p F (r) [p B (r)] . 
We assume that the fermions are not paired into a superfluid, which may require the system to have a temperature 
above the critical temperature T c for mediated interaction induced fermion pairing 0,0|- I n practice, however, the 
T c for such fermion pairing [lj| [l7| is much lower than any other relevant energy scale so that there is a temperature 
regime in which for all practical purposes we can treat the unpaired fermion-BEC mixture as a zero-temperature 
system with a normal fermion fluid. Whether the p-wave pairing can be achieved may also depend on other factors: 
for instance, the inevitable loss processes that give a finite lifetime to the cold atom trap systems also heat up the 
system fl8j . This is a process that will compete with the formation of any ultralow-temperature phase of matter. 

As mentioned before, the mean-field description commonly used in cold atom physics neglects correlations and 
fluctuations in the expression of all expectation values. Accordingly, the interaction energy of this system, -E; n t, is 
approximated as 

£int = Xbf J d\ p° F (r) p B (r) +^ j d\ p% (r) p% (r) , (1) 

where Xbf and Xbb are the interaction strengths with Xbf = ink 2 (nip 1 + m~ B l )aBF and Xbb = (47rfi 2 /tob) <Zbb, 
where asF and asB are the fermion-boson and boson-boson scattering lengths, respectively. We assume that the 
effective boson-boson interactions are repulsive, i.e., asB > 0, in order to have the system be mechanically stable. 

Molmer [191 ] pointed out that the trapped boson-fermion quantum gas mixtures can take on distinct spatial ar- 
rangements at low temperatures. When the fermion and boson particles mutually repel each other, Xbf > 0, the 
mixtures can take on spatial configurations in which one region of space is occupied by only fermion particles, an 
arrangement akin to that of the phase-separated state in condensed 3 He — 4 He mixtures. On the other hand, when 
the fcrmions and bosons attract each other, Xbf < 0, the mutual attraction causes a density increase in the regions 
of fermion-boson overlap. If the attraction is sufficiently strong, the density increase leads to a collapsing instability 



[201 ] that was observed experimentally by Inguscio and co-workers [21 1 . 

In the mixtures of repulsive fermion-boson interaction, Viverit et al. [B[ calculated the zero-temperature phase 
diagram in the mean-field approximation. Specifically, they discovered that the fermion-boson densities for an initially 
homogeneous mixture would take one of the following phases: (a) remain homogeneously mixed (phase I), (b) spatially 
separate into spatial regions of all fermion-gas and spatial regions containing a fermion-BEC mixture (phase II), and (c) 
separate into fermion-only and BEC-only regions (phase III) . They have determined the density regime boundaries of 
the phase diagram by equating the mean-field pressures as well as the mean-field chemical potentials of those particles 
that reside in both kinds of spatial regions. The pressures and chemical potentials follow from the zero-temperature 
"equation of state" : the ground-state energy of a homogeneous mixture. 

The mean- field equation of state of a homogeneous mixture of Np fermions and Nb bosons confined to a macroscopic 
volume Q, corresponding to average fermion and boson densities p a F = Np/fl and p° B = Nb/Q, respectively, is then 
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where kp is the Fermi momentum kp = (&tt 2 p F ^ . Then, the fermion and boson chemical potentials, ps = dE/dNs 
and p F = dE/dNp, follow from Eq. {2j), 



Mb = XbbPb + ^bfPf- 



(3) 
(4) 



We obtain the same chemical potential equations by minimizing the free-energy function F = E — pbNb — PfN f 
with respect to Np and Nb- Note that the homogeneous mixture represents the physical ground state only if Eqs. 
(J3j> and (J4J yield the global minimum of F. This global minimum condition requires that the second derivatives of F 
(or E) with respect to Nb and Np be positive, i.e., dps/dp B > and dp F /dp F > 0. Physically, these mathematical 
conditions imply that the compressibility of the fermion and boson systems should be positive, as the isothermal 
compressibility, kj, is inversely proportional to the chemical potential derivative, 
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where Pj denotes the partial pressure experienced by gas j (j = B or j = F), and where the subscript T indicates 
that the pressure derivative should be taken at constant temperature (zero temperature in this case). In addition, 
the condition of a global minimum also implies that the extremum cannot be a saddlcpoint, which also requires 
(SpB I 'dp F ) < {dpB/dp° B ) {dp F /dp F ). Viverit et al. 0] noted that when the fermion density exceeds a critical 
value, pF.crit, where 
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the extremum is a saddlepoint and the system can lower its energy by rolling down the saddlepoint in different Nb /Np 
directions in different spatial regions. The "rolling down the saddlepoint" suggests a mechanism by which the mixture 
can phase-separate into regions of roughly equal size, i.e., undergo spinodal decomposition. 

On the phase diagram in the boson fermion density plane, the line of spinodal decomposition, pp = pF.crit, is 
located in the region with fermion densities above the phase I /II boundary 0]. Thus, one might wonder how the 
mixture phase-separates if the densities of the initial homogeneous mixture were located in between the I /II boundary 
line and the line of spinodal decomposition according to the mean-field description. We gain insight into this dynamics 
by plotting the free energy density T ((/>, pp, ps) = F/Q, as a function of the superfluid BEC order parameter, which, 
if all bosons are Bose-condensed, is related to the boson density as |<i>| 2 = p B . Actually, we find it useful to scale 
both the BEC order parameter as well as the free energy density. We introduce |$°| with |$°| 2 = pb/^bb, and set 
(j) = $/|<i> |. We also define the scaled free energy by /, where / = FXbb/Pb- Then, substituting the fermion density 
in terms of the fermion and boson chemical potentials, we obtain 
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where the dimensionless parameters, olmf and Pmf, depend on the interaction parameters as well as on the chemical 
potentials, 
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and we have introduced an effective Fermi momentum, k F (fip) = y/ (2m p/ ft, 2 ) pp. In the density regime of interest 
(enclosed by the I/II phase boundary and the line of spinodal decomposition), the free energy Eq. ([7]) generally 
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FIG. 1: Scaled free energy / as a function of <j> [c.f.Eq. 1(7)]. Dotted line shows the stable system (one global minimum) and 
solid line shows the system that can undergo first-order phase transition (one local minimum and one global minimum at the 
origin) . 



exhibits <f> dependency as shown in Fig. Q] This shape of the Landau free energy is associated with a first-order 
phase transition: the system can reduce its free energy by climbing or tunneling the free-energy barrier that separates 
two local minima. In ordinary finite- temperature) first-order phase transitions, thermal activation locally "pushes" 
the system over the barrier nucleating clumps of matter in the new phase. If the temperature is too low to initiate 
thermal activation, the system can penetrate the barrier and nucleate by quantum- mechanical tunneling 10]. The 
latter process, however, involves many-body tunneling, which the tunneling rate generally scales exponentially with 
the number of boson particles to be formed in each clump. A recent study found that the tunneling rate exceeds the 
experimental lifetime of the system except very close to the line of spinodal decomposition or for very strong fermion- 
boson interactions. Otherwise, we can expect that at sufficiently low temperatures and below the line of spinodal 
decomposition, the fermion-boson mixture can remain in its homogeneous state for the duration of its experimental 
lifetime even if that homogeneous state represents a metastable equilibrium. 

We are also interested in what the mean-field approximation predicts for the phase-separation dynamics in the 
spinodal decomposition region, p F > pp CI - lt . We can understand the linear onset of the instability dynamics from 
the study of the collective oscillations 2| 2j|. Our studies 24 1 revealed that two collective modes can be excited: a 



zero-sound mode with its sound velocity modified by boson-mediated fermion-fermion interactions, and a BEC sound 
mode modified by the fermion mediated boson-boson interactions. As the system approaches the line of spinodal 
decomposition, the fermion-mediated interactions reduce the BEC sound velocity. This velocity vanishes at the line 
and becomes imaginary above it, signaling the exponential growth of the long-wavelength modes. We have estimated 
the time scale on which the corresponding fluctuations (which drive the onset of the phase dynamics) grow [25| . 
This analysis suggests that the fermion-mediated interactions make the BEC unstable. The vanishing of the sound 
velocity at the line of decomposition also suggests the divergence of the boson compressibility, as the compressibility 
is inversely proportional to the square of the sound velocity. However, the mean-field equation of state does not 
exhibit this behavior. The mean field also fails to include polaron effects: the mass of a fermion particle immersed in 
a phonon-fluid can be signincantlyaltered by its interactions with that fluid [2tjj ]. The importance of polaron effects 
was illustrated by a recent study Q that showed a single impurity (a distinguishable atom immersed in a BEC) can 
self-localize into a small polaron-like state. 

The rest of the paper deals with these problems by showing that the correlation physics which we will include 
in the calculation of the equation of state does lead to a divergence of the boson compressibility. The perturbation 
treatment that we develop in this paper includes polaron effects. 



III. LINEAR-RESPONSE DESCRIPTION OF FERMION-MEDIATED BOSON-BOSON 

INTERACTIONS 



In this section, we derive the fermion-mediated interaction effects up to lowest order by applying linear-response 
theory to a mean-field state. Writing the energy response of the homogeneous many-body system as the small- 
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amplitude variation of the boson and fermion static densities, we obtain 

PB (r) = p% + Sp B (r) , (10) 
p F (r) = p° + Sp F (r) . (11) 

We also determine the corresponding variation in the boson-boson interaction energy, SEbb, int, consistent with Eq. 
(fTj). To lowest order in the small-amplitude variation, the variation in interaction energy takes the form 

5E BBM t = Jd 3 rJ d 3 r' p° B (r)v BB , tot (r - r') 6p B (r') , (12) 

and from this equation we will determine the "total" interparticle interaction potential, v t ot (r— r'), that includes 
the fcrmion-mediated interaction potential, VBB,fmed(r — r'), in addition to the usual boson-boson contact interaction 
potential. 

The boson density variation causes a modification of the overall mean-field potential experienced by the fermion 
particles, 5Vf (r), where 

SV F (r) = XbfSpb (t) . (13) 

In linear-response theory, the fermion density variation 5pp (r) caused by the static fermion potential SVf (r) is 
determined by the static response function of the noninteracting fermions, Xf' 

S PF (r)= J d s r'x°F (v-t')SVf(v') 

= X BF J d 3 r'x F {v-v')8p B {v'). (14) 

The above response function is the momentum Fourier transform of the dynamic response function Xf (P) u ) m the 
limit of vanishing frequency, 

X° F (r - r') = J d 3 p x°f (P; w = 0) exp [ip • (r - r')] . (15) 
In general, the dynamic density-density response function of a system confined to a macroscopic volume Q is 



Xf (p,w) = -^IfMIpf (p)|0) F 



1 1 



fiu> - (E int - E Q ) + it] hw + (E int - E Q ) + ir) 



(16) 



where p F (p) denotes the Fourier transformed fermion density operator, Eq represents the energy of the unperturbed 
ground state |0), and the summation runs over all intermediate excited fermion states |int) of excitation energy 
E- m t ~ Eq. At the end of the calculation, we take the limit 77 — > while approaching zero from the positive side rj > 0. 
In the static limit of a noninteracting fermion system for which the unperturbed ground state corresponds to a filled 
Fermi sphere and for which the intermediate states are particle-hole excitations, x°f (Pi w = 0) yields 

X F (P; - = 0) = -1 £ ^ _ F J , (17) 

" |int) £ k+P £ k 

where e£ is the excitation energy of a single fermion of momentum k, and where denotes the zero-temperature 
occupation number of the single-particle k momentum state. Then, the interaction energy caused by Sps and SpF is 
equal to 

SE int = J d 3 r p° B (r) XbbSpb (r) 
+ J d 3 r p% (r) XbfSpf (r) 

+ j d 3 r p% (r) XbfSpb (r) . (18) 
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Inserting Eq. (fT4|) into the mean- field energy of Eq. (| 18[) . we obtain the corresponding modification of the interaction 
energy, 

SE int = J d 3 r p° B (r) XbbSpb (r) 
+ j d 3 r p F (r) X BF 8p B (r) 

+ | d 3 r | dVp° s (r) A| fXf (r - r') <5p B (r') . (19) 

A comparison of the first and third terms of the right-hand side of Eq. (fl9|) to Eq. (fT2|) indicates that we can associate 

VBB,tot (r - r') = A BB <5 (r - r') + UBS,fmod (r - r') (20) 

with the linear-response fermion-mediated interaction equal to 

vbb Jmcd (r - r') = \ 2 bfXf (r - r') . (21) 

To determine the fermion-mediated boson-boson energy contribution to the equation of state of the fcrmion-BEC 
mixtures, we integrate the fermion-mediated boson-boson interaction potential over the boson density, 

T-i f j3 f j3 / / \ v BB,fmed (r — r') Q 

E B B,fmed ~ dTr tfV p% (r) : '-p B (r') 



= %&Cp-0; W = 0)^f , (22) 

a result that we will obtain rigorously from a perturbation treatment. 

While we can evaluate Xf (p ^ 0;cj = 0) directly from Eq. (fl7|) , we can also derive it from the Thomas- Fermi 
approximation as the static (ui — » 0), long-wavelength (|k| <C kp) limit coincides with the regime of validity of the 
Thomas-Fermi description. In this approximation, we introduce a local Fermi momentum kp (r) whose local Fermi 
energy equals the difference of the system's chemical potential and effective Fermion potential energy, 

%2 [k / (r)] =Hf- XbfP° b - ^bfSpb (r) . (23) 
2rriF 

Relative to the equilibrium momentum k F , and its corresponding Fermi energy h 2 (kp) jimp = pp — \bfP%, the 
Fermi-momentum variation Skp (r), where kF = k F + Sk F , becomes 

Skp(r)^-^X B pSp B (r). (24) 
Using pf = kp/ (67r 2 ), the corresponding linear variation of the fermion density 6pF (r) is 

5p F {v)^^ r 8kp{v)=- r ^\BF&PB(r), (25) 
from which we can extract the long wavelength static density-density response function 

(26) 

With this expression and from Eq. (|22[) , we obtain the fermion mediated boson-boson interaction energy contribution 
to the mixture's equation of state, 



N 



1 _i_ a BFk°p 



Eb B.fmcd — ~^BF~^ 1 ^ ' (2?) 

2S 2 \ TUB / 7T 

which gives the fermion-mediated boson-boson interaction energy contribution that is consistent with linear response. 

We note that the fermion-mediated energy in Eq. (|27p not only arises in response to a boson density variation but 
also is an integral part of the many-body energy, as we will show in this section. 
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IV. FERMION MEDIATED BOSON-BOSON INTERACTIONS IN SECOND-ORDER PERTURBATION 

In this section, we determine the fermion-mediated boson-boson interactions using a perturbation approach. 

Two assumptions were made. One is that the mixture is homogeneous, i.e., the average densities p B and p F are 
position independent, an d the ot her is that the boson-boson scattering length, clbb, is sufficiently small to ensure 
that the gas parameter \J p B a BB < 1 so that the BEC is dilute and well-described by the Bogoliubov approximation. 
The zeroth order corresponds to the limit Abb —> while Abb remains constant. Thus, the zeroth-order ground state 
is a product state |0) = |0) s <g> |0) F , where |0) B denotes the weakly interacting BEC ground state, and |0) F denotes 
the ground state of an ideal gas of single spin component indistinguishable fcrmions, corresponding to a filled Fermi 
sphere of radius kp. Similarly, the excited states of the zeroth-order Hamiltonian are product states of the zeroth 
order fermion ground-state and the BEC ground state (or alternatively particle-hole states and boson quasiparticle 
states). 



A. The perturbation treatment 

The perturbation Hamiltonian, H p , is the effective interaction potential that describes the short-range boson- fermion 
interactions, 

4 = ^£p^ k , (28) 

k 

where p k and are the momentum Fourier transform of the boson and fermion density operators, respectively. 
Using the Bogoliubov approximation for weakly interacting bosons, we write 

Pk « N B (c B k + CB,- k ) , (29) 

where c 1 * (c) is the boson particle creation (annihilation) operator that is Bogoliubov transformed into quasi-particle 
(phonon) operators, and 77, where t]\0)b = 0. The Bogoliubov transformation also yields the energy cost of exciting 

a single quantum of the collective BEC oscillation given by the Bogoliubov dispersion hu>j? = hkc-\J 1 + (fc£) 2 , where 

c denotes the sound velocity of the unperturbed BEC, c = -\/Abb/5 b /tob and £ represents its coherence length 

£ = (lQirpgCi b b) • A direct application of the Bogoliubov transformation also gives the excitation density matrix 
element, 



B^\m B = f^f^VN-B. (30) 

The many-body energy can be calculated by expanding the term in the fermion-boson interaction strength, 

AE = E - E Q = AE 1 + AE 2 H , (31) 

where the jth-order contribution, AEj, varies as <~ (Xbf) 1 as the boson-fermion scattering length is modified. 
The first-order energy contribution gives the mean-field fermion-boson interaction energy, 

AE 1 = B (0\ F (0\H P \0) F \0) B = Abf^^. (32) 

Hence, since we can neglect the depletion contribution to the many-body ground state of the zeroth-order BEC (which 
gives a relative error of order ~ \/ Pb bb)' ^ ne zeroth-order and first order energy terms give an equation of state 
that is identical to the mean-field equation of state. Thus, the perturbation terms of higher order add correlation 
terms to the mean-field equation of state. 

The second order contribution can be written as the sum over all intermediate states, |int), 

AB2 _ E K|M (33) 
f— : E int - E a 

|int> 

where Ei nt denotes the corresponding (zeroth-order) excitation energy of the intermediate state. Here we distinguish 
three kinds of intermediate states: The first category consists of products of the fermion ground state and quasiparticle 
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boson excited states. The second category comprises products of the BEC ground state with particle-hole fermion 
excited states. The third category consists of products of excited fermion and excited boson states. It is the second 
category that yields the fermion-mediated boson-boson interaction energy, AE B ' B fmed , 



AE 



(2) 

SB, fmed 



^FB 

7T 



E! 

k,p 



*(o\pZ\oh 



(34) 



equal to Ebb, fmed m Eq. (122p . In the limit where the system size becomes infinite, the matrix element tends to 
b(0\p^\0)b — * NsSk, implying that the corresponding process involves zero momentum transfer. Perhaps because 
of the zero momentum nature of the momentum transfer processes, this contribution has been left out in other 
papers [HI, 0, EH- Here, we simply mention that this term needs to be included and we show explicitly how the 
long-wavelength limit corresponds to the infinite size limit of the fermion system. 

Consider a large (i.e., a linear size that significantly exceeds the BEC coherence length) but finite-size BEC system 
immersed in a homogeneous infinite Fermi sea. Then, let the BEC size approach that of the homogeneous fermion 
system. A finite BEC has a density expectation value f(0\p^\0)b that is a smooth function of k. A Fourier transform 
of the average spatial distribution gives 



'|Pk|0>B = J d 3 r exp^k-r) p%(v). 



(35) 



Inserting the expression of Eq. ([35| into Eq. (f34|) , we obtain 

d 3 k 



AE 



(2) 

BB, fmcd 



<fV / dV p B {v) 



\ 2 



1 ( 1 — n„ , i, 

r., / Ml 1 P V P+k 

exp [ik ■ (r - r )] - 



p P+k k 



= Jd 3 r J d 3 r> p° B (r) 



I ( n ^BB, fmed ( r r ') 



PB (r'). 



(36) 



This expression is independent of the BEC size and should clearly be included in the many-body energy regardless 
of the size of the BEC. The infinite BEC limit simplifies the expression by allowing a straightforward substitution to 
center-of-mass and relative coordinates, equivalent to the long-wavelength limit 



E 



(0) 

BB, fined 



— — hm v 

20, k^o 



BB,fmcd\^h 



(37) 



where 



(o) 

7 SB,fmcd 



(k) is the Fourier transform of ws_B,f m ed(x) in Eq. ([22 



It is interesting to note that by adding the fermion-mediated boson-boson energy to the mean-field ground-state 
energy Eq. ([2]), one obtains an equation of state that exhibits an instability at the same fermion density pF,crit as 
the mean-field equation of state. However, the instability is not a saddlepoint instability but one of diverging boson 
compressibility. The boson chemical potential derivative is equal to 



dpi 
dpi 



Abb 



^bb 



1 - 
1 - 



BF 



^BB 



1TIf\ kpd fb 



rriF 
m B 



TUB/ 7T 

ms\ Qfb flBffcf 
WW a B B 27r 



(38) 



which vanishes as pp — » PF.crit- Since the boson compressibility kb — {p°b) ®Pb /dps , the compressibility 
diverges at that same fermion density. That large values of the compressibility can make a significant impact on 
the physical behavior of systems is illustrated by the phenomenon of critical opalescence-the sudden increase in light 
scattering that can turn the appearance of systems milky white near the critical point -which is caused by the increase 
of density fluctuations that accompanies the large value of the compressibility. In trapped systems, we expect the 
increase of compressibility to result in higher densities and smaller BEC size. The strikingly different density profile 
could be an easily detectable "beyond-mean-field" effect. 



9 



boson 



fermion 



I 

I 
I 
I 




A 



P 



particle 



hole 



X interaction 



FIG. 2: Symbols and directions used for the diagrams in this paper. 



V. LINKED CLUSTER EXPANSION OF THE GROUND-STATE ENERGY OF FERMION-BOSON 

MIXTURES 

In this section, we develop a linked cluster expansion that complements and extends the above second-order result. 
The starting point of the expansion is Goldstone's theorem [27(, which states that the difference of the exact many- 
body energy, E, and the zeroth-order energy, Eq, can be written as an infinite sum of zeroth-order ground-state 
expectation values, 



where the subscript L indicates that only "linked" diagrams are included (unlinked diagrams factor out to cancel the 
normalization factor in the denominator). The expectation values of the perturbation Hamiltonian in Eq. (|39j) are 
obtained by inserting the completeness relation as a sum over zeroth-order eigenstates in between each perturbation 



operator, H p , and propagator operator, I H — E ) . Since the bra and ket state matrix elements of the propagator 



are eigenstates of Hq, the matrix elements are diagonal and give rise to unperturbed energy denominators. The 
interaction elements in the numerators (the matrix elements of H p ) are nondiagonal but are easily calculated. We 
envision the order in which the operator matrix elements occur to correspond to an effective time ordering (right to 
left corresponds to increasing "time"). We represent each term of the series of interaction matrix elements and energy 
denominators (propagators) as a diagram where the effective time parameter runs upwards (see Fig. [2]). 

In Fig. [31 a "cross", x, denotes a vertex that gives rise to an interaction matrix element. The crosses are connected 
by fermion and/or boson lines that represent the intermediate states that propagate from one interaction matrix 
element to another and give rise to an energy denominator. The excited fermion states are represented by solid lines 
with arrows: upward arrows indicate "particle" states and downward arrows "hole" states. The boson quasiparticle 
states are represented by dashed lines. Although we do not show the phonon time arrows explicitly in our diagrams, 
we treat the dashed phonon lines as upward-pointing (boson excitations do not have holes associated with them) when 
imposing momentum conservation at each vertex (see below). Then the propagator energy denominators are sums 
over particle-hole fermion excitation energies and boson quasiparticle energies. 

Since each application of the H p operator of Eq. (|28p creates and annihilates a fermion particle and either creates 
or annihilates a boson quasiparticle (or gives rise to the diagonal ground-state boson density matrix element), each 
x vertex has to have a fermion line that arrives and one that leaves, and each vertex has to have a boson line that 
either arrives or leaves. The exception to vertices that have a boson line either leaving or arriving are those vertices in 
which a boson line curves back to make a full circle-a 'loop'-which corresponds to the diagonal ground-state density 
expectation value b (0 \ p B \ 0)„ that occurs in the long- wavelength limit arising in the second-order evaluation of the 
fermion-mediated interaction discussed above. Similarly, a full line that loops back upon the x vertex that it left 
from indicates a fermion ground-state matrix element, p (0 \p^\ 0) F — * S^Np. 

The rules for translating each diagram to a numerical factor are straightforward, so we briefly describe the procedure 
here. In our case of a homogeneous system, the particle and quasiparticle excitations correspond to good momentum 
quantum numbers and we assign momentum labels Pi,P2, •■■ to the fermion lines and ki,k2,... to the boson lines, 
ensuring conservation of momentum at each x. Next, we divide the vertical 'time' axis into separate time intervals 




(39) 
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(a) 




(b) 






(d) 



FIG. 3: First- and second-order Goldstone diagrams for a single spin component fermion-boson BEC mixture, (a) Contact 
interaction, (b) fermion-mediated boson-boson interaction, (c) boson- mediated fermion-fermion interaction (direct term), and 
(d) boson-mediated fermion-fermion interaction (exchange and polaron terms). 

between the different interaction x points and determine the propagator energy denominators to multiply. Converting 
the remaining sums over boson kj and fermion momenta to dimcnsionlcss integrals introduces dimcnsionlcss 
expansion parameters a F and as that quantify the relative importance of fermion and boson excitation induced 
correlations. Then, we assign a single factor \\bf\/& of the H p operator to each of the momentum summations that 
remain after enforcing conservation of momentum, and multiply and divide by n B for the boson sums and by (i F for 
the fermion momentum sums. 

After taking the infinite system limit SI -1 ^ k — > (2ir) 3 J d 3 k, we introduce dimensionless momentum variables 

k'j = £kj for bosons, where £ = (\&irp B a BB ) is the BEC coherence length, and p- = pi/k F , where k F is the Fermi 
momentum, and h 2 kp/ (2mp) = fi F . Then integrate over fermion momenta, and write the zeroth-order fermion and 
boson chemical potentials as n F and /j, b , respectively. 
For boson momentum integrals we obtain 



1 IA 



BF 



Ms 



f\^BF\r 3 \ i 



a B ± fd 3 k', 



d 3 k' 



(40) 



where as denotes the boson expansion parameter 



olb 



1 + 2 \^bf\ 

rap J it £ 



(41) 



and for the fermion momentum sum, 



1 |A 



BF 



n 



Vf 



(42) 



where 



a _ 1 \^BF\kp _ A | m F \ |a_F_B|fc_F 

F fl° F 47T 2 V m B ) 7T ' 



(43) 



which plays the role of fermion expansion parameter. 
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The energy denominators obtained in this procedure contain fermion and boson excitation energies. So, first, we 
make the denominators positive valued, introducing a —1 factor for each propagator, and verify that particle excitation 
energies appear with a positive sign and hole excitations with a negative sign. Then, we make the propagators 
dimcnsionlcss by bringing out an energy scaling factor, ^tP F or /j, B , 



1 1 1 



e k 2fi% k'Vl 



(44) 



For those propagators that contain both boson and fermion excitation energies, we have to choose which energy 
scale to use. Scaling by the boson chemical potential, we obtain the following expression for the propagator that 
corresponds to the excited state with one BEC phonon excited of momentum — k, a fermion particle excitation of 
momentum p + k, and a fermion hole excitation of momentum p, 

1 _ 1 1 

eB k + £ p+k - 4 ~ 2 ^B k'VTTk^ + (m B /m F ) (k' • k' + 2ap' • k') ' (45) 

where k' = £k, p' = p/kp, and a = £,kp. Hence, the mixed fermion-boson propagators give rise to dimensionless 
integrands that depend on the mass ratio (mp/mB) as well as on the a parameter that compares the BEC and 
fermion lengthy scales: a = 2 (ms/nip) (ap/as)- 

We now include the contributions that stem from the interaction matrix elements. Each fermion loop contributes a 
factor Np, each boson loop contributes Nb, and each boson line that connects two vertices and carries a momentum 
kj gives rise to a factor 




N B 3 'I ' = N B \ , 3 _ ) . (46) 



Each fermion hole line of momentum p^ gives a factor . , where n F denotes the zero-temperature Fermi-Dirac 
distribution function. In scaled momentum units, 

<^n£,=0(l-Ml), (47) 

where represents the usual Heaviside function: (x) = 1 if x > and 9 (x) = if x < 0. Likewise, each upward 
fermion line of momentum pk gives the factor 

l-<-l-<, =0(|p' fe |-l), (48) 

associated with a particle fermion excitation. 

Finally, we have to specify how to take the infinite system (long-wavelength) limit associated with those diagrams 
that have loops. From the calculation of the fermion-mediated boson-boson interaction in the preceding section, we 
deduce the following rule: We assign a fictitious momentum p^ with each fermion loop and a momentum k^, for the 
boson loops, conserve momentum at each vertex with this additional momentum, and then take the long-wavelength 
limits limfci — > 0, limp^ — > at the end of the calculation. 

In the next two sections, we use the procedure described above for calculating the equation of state and then analyze 
the results. 



VI. LOW-ORDER DIAGRAMMATIC ANALYSIS 



In this section, we briefly describe the first- and second-order contributions to the many-particle ground-state energy 
of a homogeneous fermion-BEC mixture. Below, we describe the numerical contributions of each diagram shown in 
Figs. fo)-[3^(i) . The derivation of the full expressions is beyond the scope of this paper and will be reported elsewhere. 
The analysis of this section serves to reveal trends, such as the vanishing of specific diagrams in the low-density BEC 
limit cxb — » in which BEC fluctuation-induced correlations can be neglected. 

Much of the challenge in calculating the numerical contributions goes into making judicious choices in dealing with 
the different ways that the rules can be applied. The labeling of the fermion and boson lines with momenta can 
be carried out in different ways. For instance, in calculating the contribution of diagram (c), "the bubble diagram 
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with intersecting phonon line", one can either assign a boson momentum k to the phonon line and a momentum 
p — k to the fermion particle line or one can assign a momentum p' to the fermion particle line so that the phonon 
line acquires a momentum p — p'. The first choice results in an integral over a fermion momentum and an integral 
over a boson momentum, whereas the second choice results in a double fermion momentum integral. Also, one can 
scale the mixed fermion-boson energy denominators in the propagators either by the zeroth-order fermion or by the 
zeroth-order boson chemical potential. Neither of these choices affects the final result, but they can obscure common 
factors and common limits of different diagrams. For instance, diagrams (b) and (c) give contributions that are equal 
in magnitude and opposite in sign in the limit a — > 0. 

We now list and describe the contributions of diagrams (b) — (d). The diagram (b) describes fermion-mediated 
boson-boson interactions 



AE (b) = AE 



(2) 

BB, fined 



i 



rap 
m B 



CLBFkl 



(49) 



obtained earlier in Secs. lIIII and llVI In a Hartree-Fock-like analysis of the boson-mediated fermion-fermion interactions, 
diagrams (c) and (d) play a special role. Diagram (c) corresponds to the direct part of the boson-mediated interaction, 



AE {c) - AEp%, AD 



N F 

2n x 



FB 



FB 



BB 



(50) 



and part of diagram (d) gives the exchange interaction contribution. Specifically, diagram (d) corresponds to the sum 
of the second-order polaron energy shift experienced by the fermions caused by their interaction with the surrounding 
BEC phonon fluid and the exchange part of the boson-mediated fermion-fermion interaction, 



AE, 



(<*) 



AE 



(2) 

_FF,bmcd,X 



AE 



(2) 



od,P- 



(51) 



In a static approximation of the boson-mediated interaction, the fermions attract each other via an attractive Yukawa 
potential of range £. While the static approximation is not always valid, the effective interaction does have an effective 
range of order £ so that the a parameter quantifies the ratio of the mediated interaction range to the average fermion- 
fermion distance. The exchange interaction then depends on the a parameter and the mass ratio (which affects the 
region in which the static approximation is valid), 



AE 



(2) 

FF.bmcd,X 



N F 

2n x 



FB 



FB 



BB 



ex a; 



m,F 
m B 



(52) 



The dimensionless exchange function ex(ct; mp jmB) introduced in the above equation satisfies lima — > = 1, which 
ensures that the exchange and direct contributions cancel in the limit that the interaction range is much smaller than 

(2) 

the average fermion-fermion distance, as required by the Pauli principle. AE Fmcd p denotes the modification to the 
kinetic energy of ideal fermions with their dispersion altered by the interaction with the surrounding BEC, as described 
by second-order perturbation theory. This "polaron" contribution is often well described by a zero-momentum energy 
shift and an effective mass value. For the contact interaction, the second-order polaron contribution has to be 
renormalized. The resulting energy shift experienced by a single fermion particle is of order as (XfbPb)- The 
summation over fermion occupation numbers gives a many-body contribution 



AE 



(2) 

F.mcd.P 



N F a B (XfbPb) = N B otB (XfbPf) 



(53) 



which can be rewritten in the energy units of the boson-mediated interactions using NbUb — (3/2)a 3 Np (Xpb / Xbb), 



AE 



(2) 

Fined, P 



F x 



FB 



FB 



BB 



ep a 



mp 
m B 



(54) 



where the remaining a dependence is absorbed by the dimensionless polarization function ep. Note that the scale 
of the boson-mediated fermion-fermion interaction energy per fermion particle is XfbPp, whereas the scale of the 
fermion-mediated boson-boson interactions per boson particle is 



XFBP°B a F = XfbPb 1 



mp 

m,B 



dFBkf 



(55) 



This also implies that the ratio of the boson-mediated fermion energy scale to the fermion-mediated boson interaction 
energy is [pp/ (app B )~\ (Xbf/Xbb) = (8/37t)q 2 . In Fig.[U we show the dimensionless polaron and exchange functions 
as a function of a for a specific choice of the mass ration niB/nip = 6/7 (corresponding to 7 Li - 6 Li mixtures). 
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FIG. 4: Direct Ep., exchange Ex, and polaron Ep energies are scaled by N F /2Q.\ FB (Xfb / Xbb) 
a with mass ratio ms/mF = f. 



and plotted as a function of 



Below, we consider a specific region of the BEC-fermion mixture density space. Starting with a homogeneous 
mixture of repelling particles (clbb > 0,clbf > 0), comparable fermion and boson densities, and large fermion-boson 
scattering length (i.e., clef clbb, while still maintaining apF^F <C 1), we gradually reduce the boson density. In 
reducing the number of boson particles to zero, we cannot expect the second-order description of this section to remain 
valid. A few boson particles cannot mediate interactions for a large number of fermions and give a boson-density 
independent mediated fermion-fermion interaction. Nevertheless, in the process of decreasing the BEC density, 
remains constant whereas £ increases as £ ~ (Pb) 1 ^ 2 so that as — > and a — > oo as a — 2 (ms/wif) (cxf/mb)- 
Notice that in this dilute BEC limit, the contributions of diagram (d) vanish as ep. and ex rapidly tend to zero as a 
increases above unity. 

We could have suspected the vanishing of the diagram (ef) contribution as its expression involves a summation over 
a boson momentum so that we may expect at least a factor as- However, the a dependence of the dimensionless 
integrals complicates the analysis. Nevertheless, the numerical results shown in Fig. [4] confirm the naive expectation 
of vanishing energy contribution in the dilute BEC limit. This finding indicates a general trend, namely that those 
diagrams with BEC-phonon momenta that are summed over vanish in the low-density BEC limit. We can expect the 
diagrammatic analysis to simplify considerably in the dilute BEC region. Even though we consider a low BEC density 
region, part of this region displays interesting many-body behavior as the fermion-mediated boson-boson attraction 
competes with the usual contact-interaction boson-boson repulsion. 



VII. SUMMATION OF STRETCHED, LINEAR DIAGRAMS 



The general diagrammatic analysis becomes unwieldy as the values of the coupling parameters increase. The 
resulting complexity is a common failing of diagrammatic treatments, often limiting their usefulness in describing 
strongly interacting systems. Fortunately, in the fermion-BEC mixture, a significant simplification occurs in the low- 
density BEC limit. For that system, we define the low BEC density regime as the limit in which the BEC is dilute with 
respect to the usual (contact) boson-boson interactions, y p B a BB ^ 1> anc ^ m wmcn the bosons do not significantly 
affect the fermion particle properties (altering their effective mass, for instance, or including vertex corrections in the 
description of effective fermion-boson interactions). The latter conditions are satisfied if 



1 



m B 



I + ( aFB \ a FB/£, 1 



Tflp ) \ CLBB 



(56) 



In the limit that the fermion-boson interaction is increased, (1 + mp /mp) {cifb/cibb) 1, condition Eq. ([56")) is a 
stronger requirement than as <C 1. More generally, diagrams in which phonon momenta have to be summed over 
can be discarded in this low-density BEC limit. A large class of diagrams is still relevant, however, namely those 
diagrams in which fermion "bubbles" (particle-hole pairs) are connected by phonon propagators with momenta of 
vanishing value in the long-wavelength limit implied in loops. This is the case for the diagrams that we will refer to 
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FIG. 5: Higher-order diagrams that are summed. 



as "linear" diagrams-diagrams with two loops in which one can move linearly from one loop to the other in only one 
way, crossing all the interaction vertex signs (x) while passing alternating fcrmion-bubble and phonon segments. 

An important subclass of the linear diagrams is the 'stretched' diagrams in which the loops are placed at the earliest 
and at the latest time and in which the bubble-phonon segments proceed from the earliest to the latest time in a 
time-ordered fashion, as shown in Fig. [5) The (a'), (6'), and (c') diagram series shown in that graph follow from Sec. 
fVl (a), (6), and (c) diagrams shown in Fig. [3l by including larger and larger numbers of bubble-phonon insertions in 
between the loops. Mathematically, each of the next insertions introduces an extra factor z = (cxf/A) {Xbf / Xbb) to 
the previous term. The z parameter quantifies the magnitude of fermion-mediated boson-boson interactions relative 
to the usual short-range boson-boson interactions. In the infinite series summation over stretched linear diagram z 
takes on the role of an expansion parameter, 



ap X 



BF 



= 1 



4 X B b 
m B 



m B \ clbf^f clbf 
mp J 8tt clbb 



The series portrayed in (a'), (6'), and (c') take on the form A (l + z + 



z z + 



state equal to 



(57) 
(58) 

A/ (1 — z), yielding an equation of 



E 



XbbNI /1-5z\ X BF N f X 



2n 



1 - z 



2ft X 



BB 



l-Z 



X B fN f N b fl + z 



2n 



1 - z 



-N, 



h 2 k 2 F 

2m p ' 



(59) 



where in accordance to the low BEC density limit we have omitted the depletion contribution (Xbf N f /2&J y/ pbo?bb 
to the zeroth-order energy E . The second-order analysis resulted in a mechanical instability that corresponds to a 
diverging boson compressibility at a fermion density equal to z = 4 in terms of the new expansion parameter. The 
present series of stretched linear diagram summation carried out in the low-density BEC limit further shrinks the 
mechanically stable region to z < -| as the equation of state of Eq. ([59")) implies that d^s/dpB = Xbb (1 — 5z) / (10z) . 
In Fig. [51 we show the boson density derivative of the BEC chemical potential, proportional to the inverse of the 
boson compressibility, as a function of the fermion density. The density unit on the horizontal axis is /9_F jCr it, the 
critical fermion density calculated in the mean-field and second-order analysis. In the mean-field equation of state, 
the derivative dpsjdpB does not depend on the fermion density as this description neglects the fermion-mediated 
interactions. This result is qualitatively wrong, and in the second-order analysis the derivative does depend on the 
fermion density but not on the BEC density. Its value vanishes at the fermion density p_F. cr ;t, corresponding to z = 4. 
Our low BEC density limit analysis in which we sum over the stretched linear diagrams gives the lowest curve shown 
in Fig. [51 still depending only on the fermion density and vanishing at a fermion density p^stret equal to 



PF.strot — [ g | PF,crit 



0.5p 



Rcrit, 



(60) 



roughly half of the previous value. We expect that a factor of 2 difference in density should be measurable in cold 
atom traps. 

However, the stretched diagram analysis does not provide a definitive answer even in the low-density BEC limit. 
So far, we have neglected linear diagrams. The missing Goldstone diagrams can be obtained geometrically by folding 
the stretched diagrams so as to have time intervals in which multiple excitations occur simultaneously. To judge their 
importance we calculate the lowest-order folded linear diagram contributions-the corresponding diagrams are shown 
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FIG. 6: {dps / dps) (scaled with XbbT) is plotted as a function of pp (scaled with pB,crit) for mean-field approximation (dotted 
line), second order perturbation (dashed line), and infinite summation (solid line). 

in Fig. [3] Their calculation is more involved as some propagators now contain multiple excitation energies. Their 
combined contribution is equal to 



The simultaneity of phonon and fermion-particle-hole excitations in the diagrams gives rise to a dependence on the 
velocity ration r = c/vp — mpj (2msev). Note that the dependence on fermion and boson numbers is the same as 
that of the corresponding stretched diagram and so is the order in z, as well as the sign of the contribution (the latter 
only depends on the order of the diagram) . 

We expect that the main contributions to the BEC compressibility value stem from the diagrams (&') of Fig. [5] 
and their folded versions. The lowest-order folded (b') diagrams are folded version of the second diagram of the (6') 
series-the diagram that has two boson loops, two fermion particle-hole bubbles and one phonon propagator. These 
fourth-order folded diagrams contribute negative terms to the equation of state and would tend to further lower the 
value of the fermion density of diverging boson compressibility below p_F,stret- As their stretched counterpart, they 
contribute terms that are proportional to Ng and to z 2 . Unlike the stretched linear diagram series, these contributions 
will also depend on the velocity ratio r = c/vf- Hence, we expect them to correct the boson compressibility at order 
z 2 and higher. The corrections will then give a boson compressibility that depends on not only the fermion density 
but also on the BEC density. Actually, at very low values of the boson density, this dependence may be significantly 
affected by other terms as well. Because of the folded diagram r dependence, the third-order diagrams shown in Fig. 
[5] contribute a term to dps/dps that, for low values of r, r <C 1, varies as Xfb (pb/pf) zr [1 — In (r)]. At very low 
boson densities this term would, in fact, dominate, but since (pf/pb) (c/ v f) — \/ (8/3tt) (aBBkp) (pf/pb), under 
most experimental conditions, this would happen at BEC densities that would be so low it would be challenging to 
image their profiles. 



In this paper, we have developed a perturbation treatment of the zero-temperature equation of state of single 
spin fermion-boson mixtures with particles that interact via short-range interactions. This study reveals qualitative 
failings of the mean-field equation of state. For instance, the mechanical instability associated with the spinodal 
decomposition of the phase-separation transition, which showed up as a saddle-point instability in the mean-field 
description, turned out to be a point of diverging boson compressibility. A first- and second-order calculations of 
the equation of state show that the fermion-mediated boson-boson attraction overcomes the short-range boson-boson 
repulsion at the same value for the fermion density, p_F ;C ritj at which the saddle-point instability occurred in the mean- 
field analysis. We then showed that in the low BEC density limit, the diagrammatic analysis simplifies to a study 
of the linear diagrams, of which the "stretched" diagram subclass can be summed, giving an expansion parameter 




(61) 



VIII. CONCLUSIONS 
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z that quantifies the magnitude of the long-wavelength fcrmion-mcdiated interaction relative to the regular short- 
range boson-boson interactions. The resulting equation of state yields a value of the fcrmion density at which the 
boson compressibility diverges that is slightly higher than half corresponding to a factor (4/5) 3 of pF,ait- A further 
exploration of the remaining linear diagrams shows that the "folded" diagrams contribute terms that are higher order 
in z, and we speculate that these contributions will further lower the fermion density at which the mechanical phase- 
separation instability would set in, giving a critical value of the fermion density that would now depend on the boson 
density as well. In addition to revealing the inadequacy of the mean-field description near phase separation, these 
studies have also identified the fermion-mediated interaction as the cause of the mechanical instability associated with 
the spinodal decomposition of the fcrmion-boson mixtures. In a general context, the study presented in this paper 
gives a concrete illustration of how correlations and how the competition of weak interactions can become important 
near quantum phase transitions in the presence of interactions that, measured by absolute standards, remain weak. 
The results also suggest that careful experimental measurements near the phase separation transition of fermion-boson 
mixtures can explore fundamentally interesting quantum many-body behavior. 
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